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0^ . We derive the force exerted on the background plasma by an arbitrary distribution of non interact- 

ing quasi-particles, corresponding to either collective excitations of the plasma (plasmons, phonons) 
or dressed particles (photons, neutrinos). Our approach is based on the effective Hamiltonian de- 
, scribing the quasi-classical dynamics of the individual particles in the presence of a background 

medium. We recover the usual results for the relativistic ponderomotive force of a photon gas, and 
f"*^ | we derive the force, due to weak interactions, exerted by the electron-neutrinos in a background 

medium containing electrons, positrons and neutrons with arbitrary distribution functions. Gener- 
alization to other background species and other neutrino flavors is also discussed. 
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I. INTRODUCTION 



The ponderomotive force of electromagnetic waves |]-^] is a key concept in plasma physics |t],|| and plays a central 
role in our present understanding of intense laser-plasma interactions 0. This force arises whenever a non- uniform 
oscillating electric field is present in a dielectric, and can be seen as a slow time scale effect, or the average effect, due 
, '. to some non- uniformity of the high frequency oscillations of the electric field (To) . 

In general, the derivation of the ponderomotive force is based on the analysis of the single particle dynamics of 
charged particles in the presence of the electric field or in terms of the Maxwell equations for a macroscopic media 
p0| . By averaging the motion of the charged particle over the fast time scale, corresponding to the high frequency 
oscillations of the field, the slow time scale dynamics of the individual particles can then be calculated, and the net 
effect of the electromagnetic forces acting on the particle can be reduced to the ponderomotive force. In this paper, 
we present a different approach, which allows us not only to rederive the previous results, but also to generalize, 
in a straightforward way, the concept of ponderomotive force to other physical conditions or other quasi-particles 
besides photons (e.g. dressed neutrinos), where the dynamics of a single electron is not easily described in terms of 
f^S ' the classical force in the relativistic equations of motion. 

t — ■ The starting point of our treatment is a semi-classical description of the fields interacting with the plasma. The 
fields are described by their equivalent quasi-particles, or elementary quanta: photons for the electromagnetic field, 
plasmons for the longitudinal electrostatic oscillations, phonons for the ion acoustic oscillations, and dressed neutrinos 
for the neutrino field interacting with the background medium. This is done by defining a distribution function for 



& , the particles from the field intensity. Such a description is very useful because a kinetic equation can be derived for 



the correct quasi-particles distribution. The kinetics or equations of motion for the quasi-particles in the plasma are 



described by the dispersion relation, or the effective Hamiltonian, for each one of the fields or quasi-particles. The 



effective Hamiltonian can be derived either from a classical or a quantum field theory. The key point of our formalism 
is that by knowing the effective Hamiltonian describing a single quasi-particle dynamics, we are then able to derive the 
force exerted by a gas of non interacting elementary quanta in the background plasma: the dispersion relation, or the 
effective Hamiltonian, generates the equations of motion of the quasi-particles and from conservation of momentum 
K> , (action-reaction) it generates the force of the quasi-particles in the plasma. 

Even if the term ponderomotive force was coined to describe the forces acting on a dielectric in an arbitrary non 
uniform electric field, in this paper we generalize the concept to the interaction of any non uniform field with a 
background medium. In fact, the force exerted in a background medium due to a non uniform field can also be seen 
as the pressure gradient arising due to some inhomogeneity in the quasi-particle distribution. 

We aim to achieve two objectives: to derive the ponderomotive force solely based on the quasi-particle concept, 
thus providing an easy tool to generalize and unify the concept for different physical scenarios, and to explicitly 
show the relation between the ponderomotive force and the wave action density, or the quasi-particle number density. 
In particular, we will be mainly interested in determining the ponderomotive force of an arbitrary distribution of 
neutrinos (with different flavors) in an arbitrary background medium. 

Our focus will be on the derivation of the force exerted by an arbitrary distribution of quasi-particles. In Section 
II, we present our formalism and the approximations involved in our description. The general expression for the 
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ponderomotive force is then derived and represented as a function of the effective Hamiltonian or dispersion relation 
(describing the quasi-classical dynamics of the quasi-particles) and the quasi-particles distribution function. In Section 
III, we apply our results to classical fields arising in a plasma. The ponderomotive force exerted by a gas of photons 
in a plasma is rederived in the relativistic regime, thus showing the equivalence between our new approach and the 
derivation of the ponderomotive force found in the literature. We also derive the force due to a gas of plasmons and a 
gas of phonons. In Section IV, we consider the interaction of a gas of neutrinos with a dense plasma, also containing 
other species of baryonic matter (such as neutrons, protons and positrons). The ponderomotive force exerted in the 
background medium by the neutrinos is then derived. In our derivation we assume that each one of the species in 
the medium has an arbitrary distribution function, thus generalizing previous results only valid for cold plasmas [ pT| . 
Finally, in Section V, the results of this paper are summarized. 

II. GENERAL CONSIDERATIONS 

Consider a gas of non interacting quasi-particles (QPs) in a background medium. By quasi-particles, we mean not 
only elementary collective excitations of the background medium, such as plasmons and phonons, but also dressed 
photons or dressed neutrinos. Our formalism is independent of the entities we are considering as QPs, but we assume 
that the interaction between them is negligible. This corresponds to either considering collision frequencies v qq much 
smaller than the typical time scale of the process in study(dilute gas) or simply assuming that the interaction between 
quasi-particles can be discarded. Furthermore, we consider that the dynamics of a single QP is governed by the 
effective Hamiltonian iJ c g, which is a function not only of the dynamical variables of the single QP (momentum and 
position) but is also dependent of the properties of the medium where the QPs propagate. 

Using in a direct way the results already known for the ponderomotive force of e.m. waves in a dielectric in terms of 
a gradient in the laser intensity, we could immediately propose an expression for the ponderomotive force, written as 
a function of the number of photons. This approach would not show, however, the connection between the effective 
interaction felt by a single QP and the expression of the ponderomotive force, thus preventing the generalization to 
other physical conditions. Also, the dependence of the ponderomotive force on the wave action density, or QP number 
density, would not be clearly stated for non trivial QP distribution functions. Therefore, a different path must be 
followed. By determining the free energy of the system, we can relate the changes in the free energy with the work 
performed by the force exerted by the distribution of particles, and from that derive the ponderomotive force. 

We can write the total free energy as the sum of two contributions: the free energy of the background medium 
in the absence of QPs {Fq) and the additional free energy (for a given temperature and density) resulting from the 
presence of the field elementary excitations. Therefore, the free energy is written as 

dr J „f q (r,k,t)H eS (1) 

where p is the density of the medium, T is the temperature of the medium, and f q describes the QP distribution 
function in phase-space (r,k), r describes the QP position and k is related to the QP momentum p q by k = p q /h. 
For the sake of completeness we state here the most important properties of f q : 

f dk 

n q (r,t)=g sq j ^-^/,(r, k, t) , (2) 

/dk 
— / 9 (r,k,i)ff cff , (3) 

|0(k,<)| 2 = g sq J drf q (r,k,t)H cS , (4) 

where n q is the quasi-particle number density, |-0(r, t)\ 2 is the spatial energy density, and |0(k, t)\ 2 is the spectral 
energy density, so that the second term on the right hand side of equation (|l|) represents the total free energy due to 
the presence of QPs. g sq is the QP statistical weight and it accounts for spin degeneracy of the QPs. 

We now generalize the procedure described in ref. []lo|| to electromagnetic waves. Let us assume that an isothermal 
deformation u of the background infinite medium occurs. If the distribution of quasi-particles exerts a force f over 
the medium, work 8W will be done by that force, such that 
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SW = J drf-u . (5) 

where f has the dimensions of force per unit volume. We consider that this force will act only on the electrons of the 
background medium. The ponderomotive force felt by the ions can be discarded since it is m e /rrii smaller than the 
ponderomotive force felt by the electrons (m e is the electron rest mass, and m, is the ion mass). However, inclusion 
of other background species is straightforward. SW can be related with the change in the free energy for the same 
displacement u, since SW = —5(F — F ). Going back to eq. ([!]), we easily obtain 

SF = 5F (p, T) + g sq J dr J (H eS Sf q (r, k, t) + f q (r, k, t)SH cS ) (6) 

We will assume a quasi-static distribution of quasi-particles, which means that a displacement of the background 
medium will not affect the distribution f q . Furthermore, f q is a function only of the dynamical variables. Hence, the 
second term in the integral in eq.(^|) is zero. The effective energy H e g of each QP will be affected by the isothermal 
background deformation, since H c g depends not only on the dynamical variables but also on the properties of the 
background medium. The change in H e g has two contributions: (i) S^H c f{, due to the fact that particles from the 
background medium are pushed from r — u to r, and (ii) <5(2)i? e ff, due to the change of the distribution function, or 
density, of the background medium in position r. The first contribution for SH e g is 

S(l)H oS = -U • ViJcff (7) 

For the remaining contribution, we first use the fact that the relative change in the volume element of the background 
medium is dV/V — V • u. Therefore, the change in the number density is Sn^ g = — rtbgV • u, where /ib g is the 
number density of the background particles affected by the ponderomotive force (in our case, just the electrons). In 
the same way, the change in the distribution function of the background medium particles /b g (r, p) is <5/b g (r, p, t) — 
— /b g (r, p)Vu, where p is the momentum of the background particles (we also assume that the isothermal deformation 
does not impart momentum to the medium). Therefore, S(2)H e fi can be written either as 



if H c s is a function of n-bg, or 



* (2 )ff eff = -(^) A g (r,p)V-u (9) 



T 



for a more general dependence of H c g on the distribution function of the background medium /b g . Inserting eq.(Q) 
and eq.(||), in eq.(||) we obtain 

5F - 5F (p,T) = - 9sq J dr J J?L/ 9 (r,k,t) (u • VH cS + «b g V • u) (10) 

Performing an integration by parts over the second term on the right-hand side of eq. (|l0|) , results in 

SF - 5F (p,T) = -g„ J dt J An-/,(i,l,l]VI„ + 

*./*/ (^3- V (/.CM) O t *«) <") 

Therefore, the ponderomotive force, per unit volume, acting on the medium is 

/ C^/.<'.M™--*. / (^3v(/ ? (r,k,t)(|^) T n bg ) (12) 

A similar expression can also be derived when H g depends on the distribution function of the background medium 
/b g . This will be discussed for the particular case of the ponderomotive force due to neutrinos in a plasma. For a 
linear dependence of H c g on ribg i-e. nb g (dH e g/dnb g ) = H c s, eq.(|l2|) can be further simplified to 

f (ik 

f (r, t) = - 9sq J ^-^H cS Vf q (r, k, t) (13) 

Equation dl2|), and the appropriate definition of i? ff, is the starting point of our discussion of the ponderomotive 
force due to different types of QPs propagating in a plasma. It also represents a generalization for arbitrary fields of 
the Landau and Lifshitz arguments p0| . 
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III. PONDEROMOTIVE FORCE OF PHOTONS, PLASMONS AND PHONONS 



Having determined the general expression for the ponderomotive force of a gas of non-interacting QPs, we now 
proceed by evaluating eq.(|l2|) for different classical fields or classical QPs. We first consider the ponderomotive force 
due to a distribution of photons characterized by f q = Af(r, k, t). The number of photons Af, obeying the properties 
(||J]]) can be obtained from the Wigner function of the electromagnetic field |l2j . The concept of the number of photons 
was introduced in the Plasma Physics literature in the 60's, associated with the Random Phase Approximation [fl3|| . 
We point out, however, that a proper definition of the number of photons based on the Wigner function, can describe 
any e.m. field configuration (see Appendix, for a detailed discussion of this problem). For the number of photons Af, 
the role of the effective Hamiltonian H c g is played by tko(r, k, t) = Tiuj^, where the frequency Wk is obtained from the 
dispersion relation for the electromagnetic waves propagating in the plasma. This is equivalent to assuming that each 
single photon obeys the dispersion relation for plane electromagnetic waves. Since we are considering an arbitrary 
distribution of photons, the relativistic mass correction must also be included in the dispersion relation for circularly 
polarized electromagnetic plane waves: 

^ = k 2 c 2 + (M) 

where Up e (r,t) = 47re 2 n e (r, t)/m e is the local electron plasma frequency, and n e (r,t) is the electron number density. 
7 is the relativistic mass correction factor, which is a function of the electric field intensity and the energy distribution 
function of the electrons in the background plasma. For the sake of clarity, we will not write down the explicit 
expression for 7, since it will not add any additional features to our derivation. We note that in (pT|), k, r, and t are 
independent variables. Furthermore, eq.([l4|) is valid in the limits of the geometrical optics approximation E3| i.e. 

-i>-logn e (r,t) , U » |Vlogn e (r,*)| (15) 
Zir ot Zn 

These limits establish the conditions in which this approach is valid: whenever two very different time scales are 
present, we can treat the high frequency perturbations as a gas of quasi-particles propagating in a background with 
slow density modulations. 



Using eq. (114) in eq.(|12j) and the fact that 



dH cff \ ( dhuj k \ 2 hu 2 

n b g = Tr^T w pe = ( 16 ) 



dn hg J T s \dojj e J T pe 2w k7 



and 



Vtfeff = V^ k = (17) 

2w k 7 

we obtain, after some algebra, the ponderomotive force due to a photon distribution Af: 

WM) = -f W v/j^ (18) 

\E I 2 

where g p h — 1(2) for circular (linear) polarization. For plane electromagnetic waves, Af = Sg h nhio ~ k ), with E 
the electric field amplitude, and uj (fc ) the frequency (wave number) of the electromagnetic field |T^ |. In this case, 
the ponderomotive force acting on a single electron reduces to the more familiar form: 

f*(r,t) = -^VA» (19) 

where A is the high frequency vector potential. It is straightforward to see that ponderomotive force effects can 
arise due to two different conditions: inhomogeneity in the number of photons distribution function Af and/or spatial 
dependent frequency. Present understanding of ultraintense short laser pulse propagation in plasmas is based on the 
different roles played by these two contributions |L5| . Equation ([l9]) agrees with previous derivations of the relativistic 
ponderomotive force for circularly polarized photons |l||2|| For linearly polarized plane electromagnetic waves, the 
dispersion relation (|l4j) is no longer valid; a different expression for the relativistic ponderomotive force appears ||, 
which under certain conditions reduces to ( |l9| ) . As far as we know, eq. dig) is the first derivation of the relativistic 
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ponderomotive force, using solely the number of photons concept. Also, eq.dlq) is valid for any e.m. field configuration, 
as long as TV is defined using the proper definition of Wigner function for an arbitrary e.m. field (see Appendix). 

For a gas of plasmons (electron plasma oscillations) , an expression similar to eq. (|l^) can also be derived in a similar 
manner, but now the dispersion relation is 

2 2 i o ; 2 ^bTc , 

= ^pe + 3^ (20) 

m e 

where T e is the plasma electron temperature. The ponderomotive force of the gas of plasmons, described by A/piasmon, 
acting on the plasma electrons is 

c / i\ ^pe f t A/"pi asmon 

wr,t) = -— v y jw—zr- (21) 

The ponderomotive force of the gas of plasmons can then describe the coupling of the electrostatic oscillations with 
the low frequency, long wavelength ion acoustic oscillations Jl3| . 

We now consider the ponderomotive force due to a gas of elementary collective ion excitations. Due to the limits of 
our formalism (see eq. (|l5|) ) , this gas of ion oscillations must interact with even lower frequency and longer wavelength 
plasma perturbations. In a unmagnetized electron-ion plasma, such physical condition cannot be verified. However, 
dusty plasmas support oscillations with characteristic frequencies (wavelengths) much lower(longer) than those of ion 
oscillations (ion acoustic waves or ion plasma waves) [ fl6f . The picture of a gas of ion oscillations in a dust-acoustic 
oscillation is then reasonable. The general linear dispersion relation for ion oscillations is || 

^ = fc 27 ! M + ,2 IM 1 (22) 

m l rrii 1 + 7 e «r A~ 

with Ti is the ion temperature, 7 e (7i) is the electron(ion) adiabatic index, and A e = (fcsT e /47rn e e 2 ) 1 ^ 2 is the electron 
Debye length. Denoting the number density of QP excitations corresponding to ion collective motions by A/i on qp, we 
obtain from eq.(|l2|), and using eq.(p2[), the force: 



m l A| f dk ?i-Monqp / k 2 j e k B T e 
k B T e 2 J (2tt)3 uj k \ m, 1 + X 2 e k 2 j e 



fion qP (r,i) = -r^r^fV / ( < e ^\ t ±^ ) (23) 



In the limit of jiTjJrrii <C / -f e T e /rrii, the dispersion relation (^2|) can be approximated by ui 2 ~ fc 2 7 e fcsT e /{mi(l + 
A 2 fc 2 7e)}, and eq.(|23|) reduces to 

fion qP (r,i) = "^ryV J T^aMonqp^ (24) 

Two opposite physical scenarios can now be explored. In the limit of k\ e <C 1, corresponding to ion acoustic 
oscillations, the dispersion relation is Wk — kc s , with the sound speed c s = \/^i e kBT e /mi. Thus, the ponderomotive 
force due to a distribution of ion acoustic waves, or a gas of ion acoustic phonons (iap), is 

f iap M) = ---§-V / T^^Map (25) 

2uj pi J (z7r) 



with the ion plasma frequency uj p i — uj pe ^Jrn^Jmi. On the other hand, in the limit of kX e ^> 1, the dispersion relation 
( p2| ) describes ion plasmas waves (ipw), such that uj 2 ~ lu 2 ^ In this case, the ponderomotive force verifies 



tik 3 f rfk 

fi pw (r,t) = uj pi \7 I j—-^hNi vvl - - V J 3 ft^k-Mpw (26) 



We then see that the ponderomotive force associated with ion collective excitations, eqns.(^,^6|), is significantly 
different from that associated with photons, eq.(|lg|), or plasmons, eq. (^l|) . This derivation must now be verified using 
the standard plasma physics methods. The coupling between the ion collective motions and the background medium 
can play a significant role in the (de)stabilization of dust-acoustic oscillations in a dusty plasma. Both these features 
will be explored in a future work. 
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IV. PONDEROMOTIVE FORCE OF NEUTRINOS IN DENSE PLASMAS 



We now turn to the central result of this paper: the ponderomotive force of neutrinos in a plasma. This force can 
provide the coupling mechanism responsible for the anomalous scattering of neutrinos and the consequent deposition 
in the plasma of some part of the neutrino energy released in a supernovae explosion |ll] , |l7| . This energy deposition 
plays a key role in the present understanding of supernovae explosions ]l8| . An intuitive picture of the ponderomotive 
force is immediate: if in a given region the energy density in the neutrinos is higher than in other regions, a force will 
be exerted in the background medium by the neutrinos, towards the regions of lower neutrino energy density. This 
corresponds to the physical picture of the neutrinos trying to push their way into regions of lower energy density. 

In order to derive the ponderomotive force due to the electron- neutrinos (y e ) propagating in a background medium 
we first write down the effective Hamiltonian for a single electron-neutrino in a unmagnctized background of electrons 
characterized by the electron distribution function / e (r,p e ,i) |l9|: 



dp e 



H eS = y/plci + mlc* + g se / j^V cS (r, p e , t) (27) 



with V c g given by 



V cS (r, p e , t) = gvV2G F fe{v, p e , t) 1 1 - Bf-J^ 



(28) 



where Gf is the Fermi constant of weak interaction, gv = (1/2 + 2 sin 2 9w) — 1 is the effective vector coupling 
constant in the Standard Model, 9w is the Weinberg mixing angle, \) v = Kk is the electron-neutrino momentum, 
kjy = Pi//|Pjy|, E e is the electron energy (a function of the electron momentum p e ), and m„ is the neutrino mass, 
which can be set to zero for massless neutrinos. g se = 2 is the statistical weight for the electrons, corresponding to 
spin 1/2, and appears because of spin degeneracy. Also, / e (r,p e ,i) satisfies 

n e (r,t) = g se J j!^fe(r,Pe,t) (29) 

The effective potential in Eq. (^8|) has been derived using the methods of Finite Temperature Quantum Field Theory 
p0| . Hence, our approach here is clearly a semi-classical one: we assume that the interaction of the neutrinos with 
the electrons is governed by quantum processes (included in V c r), and we take into account the Fermi statistics of the 
phase space density of the particle numbers, but the neutrino dynamics is determined by the classical Hamiltonian 
obtained using the equivalence principle, and we neglect the spins. This approximation is valid as long as changes in 
V e ff occur over length scales much longer than the neutrino de Broglie wavelength A„ = 27r/|k|, and no spin waves are 
considered. 

We first analyze the contribution of the term p e -k„/i? e to the effective Hamiltonian (|2~7j). If the neutrinos propagate 
along a precise direction, and assuming an isotropic distribution for the electrons, the integral 

r/ e (r,p e ,t)— — (30) 



averages to zero. Furthermore, if the neutrino distribution is isotropic, when integrating over the contribution of all 
the neutrinos for the ponderomotive force, a zero average is obtained once again. Therefore, this term only gives a 
contribution for both anisotropic neutrino and electron distribution functions. 

Since eq. (^?j) is also a function of the distribution function of the electrons in the background medium, eq. (^J) must 
now be employed to obtain the expression for the ponderomotive force due to the neutrinos: 



fu(r,t)=g se /t^3 /?|^3^( r > k ^)VKff- 

~ 9se Iwr I Wf^ ( /l/(r ' M) (?j?) T /e(r ' Pe " 



t) (31) 



where /„(r,k, t) is the neutrino distribution function, which also verifies the normalization condition 

f dk 

n v (r,t)=g su J —^f u (r,\s.,t) (32) 
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where n u (r, t) is the neutrino number density, and g sv — 1 is the neutrino statistical weight (neutrinos are completely 
polarized (left) particles). Since the effective potential has a linear dependence on f e (v,p e ,t), then 



df, 



f e (r,p e ,t) = V c g (33) 



e / T 



and eq.([3l|) reduces to 



fv(r,*) = -9se Jj^J 72§3^ ff V/,(r,k, 



t) (34) 



Inserting eq.(|28|) in eq.(|3J), we obtain 

v 7 ^/-, • 2/1 f dk [dp 



Ur,t) = -^f(l + 4sin 2 9 w )G F9se J ^ j T^/efr Pe, *)V/„(r, k, f) + 

y^(l+As^e w )G F g se J ^ | ^/ e (r ) p e ,t)H^V/ l/ (r,k, i ) (35) 

Using eqns.( p9| , p2[ ), the first term in eq. (|35|) can be re- written as 

/2 

f„(r,t) = -^-(1 + 4sin 2 6W)GFn e (r,i)Vn 1 ,(r,£) (36) 

Eq.(|35|) represents the force per unit volume exerted by the neutrinos over the electrons contained in the unit volume. 
The second term, which accounts for the contributions from anisotropies, is different from zero only when the electron 
distribution function and the neutrino distribution function are anisotropic. This means that even for a beamed 
neutrino distribution, the second term vanishes for an isotropic plasma. From now on we will consider an isotropic 
electron distribution function, and hence we discard the contribution of the second term in eq.(|35|). The first term in 
eq.(|3^) (or eq.(|36|)) is then responsible for the strong coupling between the neutrinos and the electrons, as suggested 
before jlTj . From eq. (^) , we can easily derive the force acting on a single electron due to the presence of the neutrino 
distribution: 

/2 

f„_ e = -^-(1 + 4sin 2 ^)G F Vn„(r,i) (37) 



We stress out that Eq.(36) is valid for any neutrino or electron distribution functions. 

The generalization of Eq. (5^ , ^ , ^) for the force exerted by neutrinos and antineutrinos over electrons, positrons or 
the neutrons is also straightforward. Changing the effective potential V e g pl[ | in order to include the weak interaction 
of the electron-neutrinos and anti electron- neutrinos with electrons(positrons), we obtain the ponderomotive force 
over a single electron(positron) 

fx 

f„p- e +(e-) =T — (l + 4sin 2 H /)Gi,VK e (r,i)-n Pe (r,i)) (38) 

where the minus(plus) sign refers to electrons(positrons), being n Pe (r,t) the anti electron-neutrino number density. 
The ponderomotive force over a single neutron can also be written as 

%v-n = —G F V (n ve (r,t) - 7ie e (r,t)) (39) 

As before, we are assuming a unmagnetized background medium, and we are discarding the contribution of the 
anisotropies of the electron/positron/neutron and neutrino/antineutrino distribution functions. Generalization of the 
ponderomotive force due to other neutrino flavors (tau-neutrinos or muon-neutrinos) is straightforward as long as the 
proper effective potentials V e s are considered 

It is now important to compare our expression of the ponderomotive force with that first introduced by Bingham 
et al in a phenomenological way, or the one derived by Hardy and Melrose (2^] using the methods of Quantum 
Plasmadynamics [ p3[ . The expression derived by Bingham et al is based on the analogy between the ponderomotive 
force due to electromagnetic waves (as derived in JlO{] ) and the ponderomotive force due to the neutrinos, and it 
is only valid in the limits of validity of the Landau and Lifshitz expression i.e. as long as the neutrino flux is 
assumed monochromatic and the energy (or frequency) of the neutrinos is assumed constant. In this particular 
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physical scenario eq.(ll) in ref. |11| is equivalent to our eg. (p6|) . A comparison with the results derived in ref. |22| ] 
shows that their results are equivalent to those presented here, in the limit of constant background electron/positron 
number density. The discrepancy is present whenever gradients of the background number density are assumed, and 
it arises from a misinterpretation of the ponderomotive force concept in reference p2~| . A correct interpretation of the 
results obtained by the Quantum Plasmadynamics formalism gives the same expression for the ponderomotive force 
of electron-neutrinos in a background of electrons and positrons as the one derived here [^| . 

It must be stressed that a long range interaction force between neutrinos and electrons was identified before in 
the context of a quantum kinetic treatment of neutrinos in a lepton plasma. An equivalent expression to eq. (|36|) is 
evident in eq.(18') of ref. [^5). The long range interaction force is associated with an effective charge the neutrinos 
acquire in a background of electrons |25|,|26| . Our results allow us to clearly identify the long range force in [^|| as the 
ponderomotive force due to the weak interaction of the neutrinos with a plasma. 



V. SUMMARY 



In this paper, we have presented a general derivation of the ponderomotive force due to quasi-particles propagating 
in a background medium. In particular, the results for the relativistic ponderomotive force of photons and for the 
ponderomotive force due to plasmons were recovered. The force of a gas of ion collective motions in the background 
plasma was also presented. We then applied the same techniques to a gas of neutrinos interacting with a background 
medium through the weak interaction force. The ponderomotive force due to an arbitrary distribution of neutrinos 
and antineutrinos interacting either with electrons, positrons or neutrons was derived. 

The ponderomotive force derived here is the force on a Lagrangian fluid element or on a single particle from the 
background medium. It then provides the proper way to build a single-particle or a self-consistent kinetic theory for 
the interaction of QPs (neutrinos, photons, plasmons, phonons) with a background medium, which can be applied as 
the foundation tool for the study of QPs driven instabilities in a plasma. 
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98DP00211 and DE-FG03-92ER40727. LOS acknowledges the financial support of PRAXIS XXI through grant 
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VI. APPENDIX 



In order to clarify the meaning and scope of the number of photons concept, we present a calculation of the number 
of photons Af for several e.m. field configurations. This will allow us to connect this concept with previous definitions 
found in the literature and to clearly point out that no Random Phase Approximation is assumed when a proper 
definition for the number of photons is employed. 

The number of photons Af obeys the general properties described by eqns.fUJi]) for the distribution function f q . 
The class of phase-space distribution functions that verify these properties for a given wave field, are usually denoted 
as Wigner functions |p7|] . The most common representation of the Wigner function for the electric field E(r,t) is 

Af{k,r,t) = - -. / dsE(r- s/2,i) • E*(r + s/2, i) exp(ik • s) (A.l) 

og p hTrnuJk(r, t) J 

where Wk(r,i) is obtained from the dispersion relation D(w,k, r,t) = 0. For the wave field described by Af, equation 
(^J) defines the wave action density, and eq.(|^) defines the field energy density. The Wigner function exactly satisfies 
a kinetic equation which reduces to a Vlasov equation in the short wavelength high frequency approximation. 

For simplicity, we will consider propagation in vacuum, thus meaning that Wk = |k|c. We now calculate Af for 
different electric fields. We first start with a plane wave monochromatic wave, E(r,i) = Eoexpi(ko • r — uJot), being 
uio the frequency of the electric field and k the wavevector. The number of photons is then simply written as 

A/- planc (k, r, t) = |E °j S(k - k ) (A.2) 

describing a monochromatic beam of photons, as expected. This is also the usual result present in the Plasma 
Physics literature since the 60's |2l| . A test of the ability of the Eq.(A.l) for the definition of the number of photons 
Af to describe the complete structure of the electric field arises whenever the electric field depicts an interference 
pattern. The most simple case corresponds to an electric field described by the superposition of two plane waves 
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E(r, t) — Eo exp z(ko • r — toot) +Ei exp i(ki >r — u>\t). The number of photons for this electric field defined by Eq.(A.l) 
is 

A/Bw(k,r,t) = 1 —— [|E | 2 ,5(k-ko) + |E 1 | 2 «5(k-k 1 ) + 

8g ph Trh\k\c 

E ■ E; cos( A fc ■ r - A u t)6 (k - k ° ± kl j (A.3) 

where Afc = ki — ko, and A w = — Wo- The presence of two photon beams, associated with the two plane waves 
is obvious. Furthermore, a third beam is present which results from the interference between the two plane waves, 
showing the characteristic slow modulation of the beat pattern. Using eq. (^) , we can immediately recover the usual 
energy density for two interfering monochromatic planes waves. It is then clear, that the number of photons, as 
described from eq.(A.l), already contains the information about any interference pattern present in the electric field. 

To clarify this point, and make a bridge to the usual definitions, we consider the superposition of several monochro- 
matic plane waves, also containing an additional random phase factor i/> k . The electric field is then written as 
E(r,i) = J2 k , Afc' expi(k' • r — uj^t + V'k')- For this electric field, eq.(A.l) verifies 

k' + k" 

x expi(A k , k » ■ r - A^t + A^//)5(k — ) (A.4) 

where Ak'k" = k' — k", A U ' U « = u>' — u>", and A^i^n = tpw — V'k"- Once again, the beat interference pattern 
observed in eq.(A.3) is also present for the terms verifying A k 'k" 7^ and A^/^" 7^ 0. Eq.(A.4) is the generalization 
of eq.(A.3) for the superposition of an arbitrary number plane waves. We stress that so far, no assumptions have been 
made regarding the properties of the phases ipk. When the phases "0k are random, a phase averaging of cq.(A.4) can 
be performed. This averaging corresponds to the well known Random Phase Approximation (RPA) pq | . Being the 
phase averaging =<> defined as the average of a statistical ensemble of systems differing from one another only in 
the phase -0k or A^/^w, it is obvious that < exp iA ,/,/,/," >= S(k' — k"), thus leading to the RPA number of photons 
distribution function 



A/R PA (k,r,i) = <A/R P (k,r,t) >= 1 ]T ^ A k ,A k „ 



8ff pfl 7r^|k|c 

x exp i(A k , k „ ■ r - A^ u „t)S(k' - k")<5(k - ^^-) 

= z Vr^y"l^k<| 2 £(k-k') (A.5) 

8g ph irh\k\c V 

This is the conventional definition of the number of photons, only valid under the limits of the Random Phase 
Approximation. No interference pattern is present, thus describing independent and non interfering photon beams. 

From this discussion, it becomes evident that a definition of the number of photons based on the Wigner function 
can rigorously describe different e.m. field configurations, including those where interference between different field 
components is important. In the limit of the Random Phase Approximation, the usual definitions are recovered. It 
may be argued that the Wigner function presents some pathologies (it is not a positive-definite function) which, at 
first sight, could prevent its use. However, we stress that the quantities with straightforward physical meaning are 
the marginals of the Wigner function, eqns.(||,||,[i]), and these possess the correct physical properties. 
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